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Abstract 

We analyze a Bose-Einstein condensate (BEC) - based free oscillation atom Michelson interfer- 
ometer in a weakly confining harmonic magnetic trap. A BEC at the center of the trap is split into 
two harmonics by a laser standing wave. The harmonics move in opposite directions with equal 
speeds and turn back under the influence of the trapping potential at their classical turning points. 
The harmonics are allowed to pass through each other and a recombination pulse is applied when 
they overlap at the end of a cycle after they return for the second time. We derive an expression 
for the contrast of the interferometric fringes and obtain the fundamental limit of performance of 
the interferometer in the parameter space. 
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I. INTRODUCTION 



Atom interferometers using cold atoms or Bose-Einstein condensates (BECs) can have 

applications in field-sensing and precision measurements [2]. Atom interferometers can be 
more versatile than the optical ones and have been used to measure acceleration [3| , rotations 
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^, and dynamic polarizability of atoms |5[. 

The first atom interferometry experiments with supersonic atomic beams were reported 
in Ref. jo, 7|. The laser cooling techniques of neutral atoms developed in the 1980s ^ 
opened up the applications of cold atoms in atom interferometry. Atom interferometry with 
cold atoms by projecting them in a vertical direction was used to measure acceleration due 
to gravity [3|. 

After the experimental realizations of BECs in dilute atomic gases in the mid-nineties 
9 111' the horizon of atom interferometry has broadened. The atoms in BECs have a very 
narrow momentum distribution and hence can be controlled and manipulated more easily 
than the thermal atoms by using light waves. Moreover, all atoms in BECs are in the 
same quantum state and hence BECs are excellent coherent sources of matter waves. The 
interference of two independent condensates was first reported in Ref. jl2j], in which two 
separate condensates were prepared in a double-well potential and allowed to interfere by 
switching off the potential and letting the condensates expand. Shin et al. [13j showed 
trapped atom interferometry with a condensate prepared in an optical single-well potential 
and then coherently split into two by deformin g th e single-well into a double-well potential. 
This, as well as several other experiments [l|,ll4-16| on BEC-based atom interferometry, 
shows that condensates are good candidates for interferometric applications. BEC-based 
atom interferometers in Michelson geometry [l, and in Mach-Zehnder geometry 18, 1^ 
were realized recently. 

The basic steps of guided atom interferometry are the following 20|]: an atomic wave 
packet is split into two in a trap or a wave guide, the split wave packets are sent down two 
different paths, and recombined at the end of the interferometric cycle. For example, in 
a single-reflection atom Michelson interferometer, a BEC in a zero momentum state ipo is 



split at time r = by a laser standing wave into two harmonics ^/'^ and ip. 



m 



The 



atoms in the harmonic absorb a photon from a laser beam with the momentum Hki and 
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re-emit into the beam with the momentum —hki (with ki being the wave number of the 
laser,) thus acquiring velocity vq = 2hki/M, where M is the atomic mass. Similarly , an 
atom in the ip- harmonic acquires velocity —Vq = —2hki/M. At time r = T/2, where T is 
the interferometric cycle time, a reflection pulse is applied to reverse the momenta of the 
harmonics. At time r = T, the two harmonics are subject to the action of a recombination 
pulse. After recombination, in general, the atoms populate all three harmonics ipo and 
The number of atoms in each harmonic depends on the relative phase acquired during the 
interferometric cycle and can be used to deduce this phase. 

Loss of contrast in a single reflection interferometer is primarily due to a coordinate- 
dependent phase acquired during the cycle. This phase is caused by the confining potential 
and the velocity mismatch due to reflection pulses. To overcome this drawback, a double 
pass interferometer with reflection pulses was proposed and implemented in Ref . [17|] . In this 
interferometer, each cloud travels in both arms of the interferometer before the clouds are 
finally recombined. Because of the symmetry in the paths followed by the two clouds, the 
coordinate dependent phase is partially canceled. But the reflection pulses still introduce a 
velocity mismatch that limits the performance of the interferometer [22|. The performance 



can further be improved in the geometry of a free oscillation interferometer [15|, |23|] that 
does not use the reflection pulses at all. In this geometry, the split clouds turn around in 
the confining potential at their classical turning points and pass through each other twice 
before they are finally recombined. 

Both double reflection and free oscillation interferometers are not suitable for measuring 
static environment like gravitational acceleration because both clouds accumulate equal 
environment-introduced phases during an interferometric cycle and hence the relative phase 
shift is zero. But they can be used for measuring environmental effects that can be controlled 
in time. For example, Deissler et al. in Ref. |5| measure the dynamic polarizability of ^^Rb 
atoms with a double pass interferometer. Burke et al. in Ref. 2^ show that a double pass 
interferometer can be used as a Sagnac interferometer to measure rotation using the Sagnac 
effect. 

In this paper, we theoretically analyze a free oscillation atom Michelson interferometer in 
the framework of a mean field approach and obtain its limit of performance in the parameter 
space. This interferometer is schematically shown in Fig. [H The solid lines are the paths 
followed by the two harmonics during an interferometric cycle and the vertical wavy lines 
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represent the splitting and recombination laser pulses. The split condensates move in a 
weakly- confining harmonic trap and are reflected from their classical turning points. They 
pass through each other, reach the maximum excursions in the opposite arms and return 
again. The harmonics are recombined when they again ove rlap at the center of the trap. This 
interferometer has been experimentally realized in Ref. 15|, |23(| and in a different (Mach- 
Zehnder) geometry in Ref. 16]. Horikoshi et al. IgI have shown that dephasing in this 



type of interferometer is suppressed. Nevertheless, there is still a fundamental limit on the 
performance of this geometry that is caused by the conflning potential and the nonlinearity 
of the condensate, which is the subject of the present paper. 

The rest of the paper is organized as follows: Section [III formulates the analytical model 
used to describe the interferometer. The equations of motion for the split condensates 
are derived and analyzed in Section IIIII The limits of performance of the free oscillation 
interferometer are discussed in Section IIVI The free oscillation interferometer is compared 
with the single and double reflection interferometers in Section |Vl Finally, conclusions are 
presented in Section IVll 

II. ANALYTICAL MODEL 



The results of this paper are obtained in the framework of a one- dimensional (ID) mean- 
fleld theory in Thomas- Fermi limit 25|. A ID model is a good approximation to the ex- 
perimental situation Q, Q-H, Q, where the BEC clouds are cigar-shaped with the 



largest dimension along the weak guiding direction of the trap and are moving along the 
same direction. 

Speciflcally, we describe the evolution of a BEC in a weakly-conflning parabolic potential 
of longitudinal frequency u by the following dimensionless Gross-Pitaevskii equation (GPE): 



^(X,r) 



(1) 



The axial coordinate x is normalized to the initial longitudinal radius Lq of the condensate: 
X = x/Lq. The dimensionless time r is given by the relation r = ut, where u is the longitu- 
dinal frequency of the weakly-confining potential. The strength of interatomic interactions 
is described by the parameter gm = 2u±asN/ {ujLq) where is the s-wave scattering length, 
is the total number of atoms in the condensate and uj±_ is the trapping angular frequency 
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in the tightly confined transverse dimensions. Finally, e = (ao/Lo)'^, where ao = ^/h/ (Moo) 
is the oscillator length along the longitudinal dimension. The wave function ip has 
normalized to 1. For details of derivation of Eq. ([T]) and its limits of applicability, see 

The initial equilibrium size of the condensate in the Thomas-Fermi approximation j25| is 
given by : 

^° - Mco^ ) ■ 

The wave function ip of the condensate after the splitting pulse is a superposition of two 
harmonics ipj^ and 
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(3) 



The wave functions il)± have been normalized to 1. 

The densities n± and the phases (j)± of the harmonics and if)- defined by the relations 
■^i = A/^exp(?0-|-), are represented as : 
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±V{XtXo] 



2\R) Q\R 



(4) 



In Eq. (jll), ±Xo(r) are the positions of the centers of the two harmonics and R{t) = L/Lq 
is their dimensionless radius. Since the splitting pulses act for a very short period of time, 
the harmonics' shape and position immediately after splitting are equal to those of the 
initial BEC at rest, i.e., Xo(0) = and -R(O) = 1. In the expression for the total phase of 
the split BECs (the second equation of the Eq. (jl])), the term (0o)± is the phase accrued 
by the harmonics from the environment. The term ±V{X =p Xq) is due to the motion of 
the two harmonics. The parameter V{t) is the normalized speed v of the harmonics, i.e., 
V = v/ (ujLo) with the initial value V{0) = Vq = Vq/ (wLq). The quadratic term proportional 
to G appears because of dispersion of the harmonics. The cubic term proportional to S 
is due to atom-atom interactions in the condensate. The quadratic and cubic phases are 
initially zero, G{0) = S{0) = and evolve with time when the harmonics start propagating. 
The parabolic form of the density of the BEC clouds n± in Eq. (jl]) implies the Thomas- Fermi 
limit, when the second derivative of the densities n± in Eq. ([T]) is neglected. 

The Gross-Pitaevskii equation ([1]), with the form of the density and phase of the two BEC 
clouds given by Eq. (jl]), has been previously investigated under various approximations both 



analytically and numerically in 15|, |22|, |26|, |27|]. In particular, the authors of Refs. [22| 



derived the set of coupled ordinary differential equations for the parameters -R(t), Xo{t) 



V{t), G{t), S{t), entering the expressions 



'or the density and the phase of the BEC clouds 



given by Eq. (gD (notations of Refs. |22|, Q are slightly different from those of the present 
paper). Validity of the analytical model has been confirmed by comparing solutions of these 
equations to the results following from direct numerical solution of the Gross-Pitaevskii 
equation ([1]). The derivation of equations in Refs. {22], Q has been based on representing 
the density and phase of the two BEC clouds (jl]) in terms of a truncated set of Legendre 



polynomials. In this paper we derive the set of equations analogous to that of Refs. |22l . [27 1 
by analyzing equations of motion for the quantum-mechanical expectation values associated 
with the parameters -R(t), Xo(r), V{t), G{t), and S{t). This approach streamlines the 
derivation and allows for a greater insight into the physics. Additionally, it has the energy 
conservation law explicitly built in the formalism, greatly assisting further analysis. Finally, 
the new approach is more readily generalizable to two- or three-dimensional case. The set of 
derived equations is used to analyze the geometry of the free oscillation atom interferometer. 

III. DYNAMICS OF THE SPLIT CONDENSATES 
A. Equations of motion 

The time evolution of the expectation value (A) = (tplAltp) = J %jj*Ailjdx of a quantum 
mechanical operator A is governed by the equation 

^{A) = tmH,Am, (5) 

where H is the Hamiltonian of the system and [H, A\ = HA — AH is the commutator. The 
Hamiltonians of the two BEC clouds after the splitting are given by the relation 

H± = '-P^ + + \giD{n^ + 2n^). (6) 

The expectation values of the coordinate X, i.e., {ip\X\ip) and momentum P, i.e., {il)\P\ip), 
evaluated with respect to the wave function are : 

{A'> = X„, 



Similarly, the expectation values of X^, P^, n+ and n_ are : 



lOPy 5P2 I75i?2 

= 4. (8) 
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(n-) = — (5gc?i + ^2), 
5/t 

where the functions di are defined by the relations (we shall need ^3 slightly later) : 
di{q) = qi\q\-lfi\q\+2)9i\q\<l), 

dM = {\q\ - lf{-Q\q? - 12|g|2 + 2\q\ + l)e{\q\ < 1), (9) 
dM = 35g|g|(|g|-l)2(2|g|2 + 4|g|-3)^(|g|<l). 

Parameter q = Xq/R in Eq. ([9]) is the relative position of the center of mass of a harmonic. 
The ^-function in Eq. (Q is equal to one if its argument is a logical true and zero if it is a 
logical false. These functions arise because of the interatomic interactions between the two 
harmonics. Therefore, they are non-zero only when the harmonics are overlapping. 

Using expectation values given by Eqs. (I7j), ([H]) and evaluating their dynamics with the 
help of Eq. ([5]), results in a set of first order differential equations describing the dynamics 
of the split condensates: 

- P' 



WRJ^ R^ 

d-i{q) 



Sr 



R 



where represents the derivative of the function A with respect to time. 

In deriving the equations of motion (Eq. fTOl) , e and Vq'^ were used as smallness parameters 
and terms of the order of and V^~^ have been neglected. This can be justified by the 
following estimate. For a BEC of ^"^Rb atoms, Vq = 11.7 mm/s . For the longitudinal angular 
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frequency = 27r x 4.1 Hz , the angular frequency in the transverse dimensions u± = 2'n'X 80 



Hz 23|, and the number of atoms in the condensate = 10'* [1], the equihbrium size of a 
condensate Lq given by Eq. is approximately 40 fim. For these parameters, the inverse 
of the dimensionless initial speed Vq of the harmonics is Vq^'^ ~ 0.09 and e ^ 0.018. 

Finally, analysis of Eq. f llOp requires evaluation of functions Di{q) = di{x)dx, which 
are integrals of the functions di{q), with respect to q: 

DM = ^g'(|g|'-5|g| + 5), 
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D2{q) = -g(|g|-i)'(g' + 3|g| + i), (11) 

DM = ^|g|^(20g^- 126g2 + 175|g| -70). 

Expressions for the functions -Dj(g) given by Eq. ( ITTj) are valid in the region |g| < 1. For 
|g| > 1, the functions Dii^q) are constant and equal to their boundary values : D\{±X) = 1/5, 
D2{±1) = 0, and Ds{±l) = -1/2. 

B. Evolution of the radius and the quadratic phase 

Time dependence of the radius R{t) and the quadratic phase G{t) can be obtained by 
solving the first two equations of Eq. ( ITOj) . The contribution from the term with d2{q) in 
the equation for can be neglected and these equations reduce to : 

^rr = -i?+^, (12) 

Integrating the Eq. (fT2|) with initial conditions -R(O) = 1 and Rr{0) = yields 



The first of Eq. ( ITOl) then gives 



G{t) = ±^/R{l-R){R^ + R-l). (14) 

The analytic solution to Eq. f[T^ can be obtained in terms of elliptic functions. It is 
important to notice that Eq. (|T3l) for R (and, thus, Eq. IHM for G) is "universal", i.e., 
independent of the trap frequencies, number of atoms in the condensate, etc., and needs 
to be solved only once. Fig. [2] shows the time evolution of R and G for a full trap period 
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obtained by solving the first two of the Eq. (fTOi) numerically. The small kinks in the plot 
of G during splitting, at recombination , and when the harmonics pass through each other, 
are due to mutual interaction between the two harmonics at overlap (the term with d2 in 
equation for Gr)- Fig. |2] shows that the neglect of this term in obtaining Eq. f lT5]l is an 
excellent approximation. It is interesting to note that the period of the oscillations of the 
radius is about 60% of the trap period. The quadratic phase G has the same period as R. In 
our analysis we shall need the values of R and G only at the time of recombination r = 27r: 
i?(27r) ^ 0.81 and G{2ti) ^ 0.27. 



C. Evolution of the cubic phase 

Evolution of the cubic phase S{q) is governed by the last of Eq. flTU]) . The cubic phase 
changes only when the harmonics overlap because, otherwise, the function d^{q) in (Eq. [H]) 
is zero. Integration of Eq. flTU]) yields 

because in the region of overlap dq/dr ^ ±Vo/R and S{0) = 0. The function D^{q) is given 
by Eq. ffTTj) . After the first separation of the harmonics, the value of 5* outside the overlap 
region is 

Sil) = (16) 

because D^{1) = — 1/2 . The difference between the values of S before and after the passage 
of the harmonics through each other around mid-cycle r T/2 is zero , since its calculation 
involves integration of the odd function d^{q) from g=ltog = — 1. Finally, near the end 
of the cycle 

S{q)-S{-l) = ^[DM-D,{-l)]. (17) 
Combining Eqs. ( !T6|) and (fT7|) . the value of S near the end of the cycle 



(18) 



in the lowest order of 

Fig. E] shows the evolution of S with time. It is zero initially and grows to a negative peak 
once the two harmonics start moving away from each other. After the harmonics completely 
separate, the value of S remains constant at its boundary value. The inset picture in Fig. [3] 
shows the evolution of 5* when the two harmonics pass through each other. 



D. Energy of the system 



The total energy E of the condensate is obtained by evaluating the expectation value of 
the Hamiltonian {iplHlip) with respect to the total wave function ip of the split condensates 
(Eq. [3]). This quantity can be represented as the sum of the three terms: 



where the d - functions are given by the Eq. iQ. The term Ek in Eq. (|T9l) is the kinetic 
energy of the system. It depends upon the speed of the clouds, the quadratic and the cubic 
phases. The kinetic energy right after splitting is proportional to Vq but, when the clouds 
evolve with time, the quadratic and the cubic phases develop and the kinetic energy has 
terms containing speed as well as these phases. The term Ep is the potential energy of the 
clouds. This energy is equal to the sum of the potential energy of the condensate in the 
trap caused by its finite size, and the potential energy due to the displacement of the center 
of mass of each cloud after splitting. Due to its finite size, a cloud has non-zero potential 
energy even when it is at the bottom of the trap. Finally, the term E^ is nonlinear energy 
due to atom-atom interactions in the condensate. 

During the interferometric time between the splitting and recombination pulse, the total 
energy of the system is conserved. 

E. Wave function at recombination 

The harmonics are recombined by using a recombination pulse at the end of the inter- 
ferometric cycle time, T. Since we are considering the case of {4>o)± = in Eq. (jl]), the 
recombination pulse, in the ideal situation, should recombine the two harmonics into one at 
rest. But, because of the spatially-dependent phases accumulated during the interferometric 



E = Ek + Ep + En, 



(19) 



with 




(1 + lOqdi + 2d2) 



lOR 
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cycle, there will be ^ three harmonics in the output ports - one at rest and two moving in 
opposite directions [22]. The wave function for the zeroth order harmonic at recombination 
is given by the expression 

V'o(e,g) = y^Ocos</>(e,g), (20) 

where the spatial relative phase across a harmonic is 

0(e,g) = Air(g)e + r(g)e^ (21) 

The strengths of the linear and cubic phases as functions of q is given by the relations 

Rt ( -r. -J. Gt , St 2 



m = 1^, (22) 

where Rt, Vt , Gt and St are evaluated at time T. In Eqs. (!20l) and ( l2Tl) . ^ = X/Rt is 
the normalized coordinate and q = Xq/Rt is the normalized position of the center of mass. 
We have neglected a small degree of incomplete overlap in the densities, but have taken it 
into account in phases. 

IV. FRINGE CONTRAST AND LIMITS OF PERFORMANCE 

The population in the zeroth order harmonic Nq is given by the expression 

No = Ntot J \M^,Q)\'d^, (23) 



where ipo is defined by Eq. (1201) . In the ideal case AK = F = 0, the contrast of the 
interference fringes V = 2NQ/Ntot — 1 should be equal to one, but the coordinate-dependent 
phase results in a decrease of the contrast. The physical reason for this decrease is explained 
by Fig. m When the coordinate-dependent phase q) becomes large, the density profile 
undergoes several oscillations between its maximum value (given by the envelope n{^)) and 
zero. In this limit, the relative population in the zero momentum state after recombination 
No/Nfot = f d^li'oiOl'^ approaches 1/2, and the contrast V = [2{NQ/Ntot) — 1] goes to zero. 
The contrast of the fringes for small values of AK and F can be expressed as 



(AKf + ^AKT +^F2 



(24) 
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Stickney et al. [27| noticed that the best contrast does not necessarily correspond to the 
complete overlap of the BEC clouds q = and sometimes can be improved by slightly 
changing the recombination time (or, equivalently, the value of q). We shall evaluate the 
limits of performance of the interferometer by minimizing the bracketed quantity on the 
right hand side of Eq. f l24p with respect to q and setting the result (somewhat arbitrarily) 
to 0.5. 

The velocity Vr of the cloud at recombination time in the expression for AK given by 
Eq. fl22|) can be obtained by equating the total energy of the clouds at recombination time 
T to their total energy right after splitting (r = 0) because the total energy of the system 
is constant. From Eq. (IT^ . the total energy of the system at time r = 

^(0) = I (v,' + , (25) 



2 

and the total energy at recombination time, T 



E{T) = ^ + + ^ + (26) 

^ ^ 2 lORl 10 IORt ^ ' 



Equating the right hand sides of Eqs. (]25l) and (1261) gives 

2Vo \^Rt 5 5i?T 5/ 
where Vt is the velocity of the harmonic right before the recombination pulse. By substi- 
tuting {Vt — Vo) from Eq. ( 1271) into the first equation of Eq. ( l22l) . one gets 

Air(g) = -i(^ + 0.27g), (28) 

where q = Xq/Rt is the relative position of the center of mass of the ip+ harmonic. To get 
AK in the form given by Eq. ( 128|) . we used Rq = 1, Rt = 0.81, and Gt ~ 0.27. Similarly, 
from the second equation of Eq. ( 122|) and Eq. (fT8|) , 

r(,) ^ (29) 

Minimizing the bracketed quantity in the right hand side of the Eq. with respect to q 
and requiring V > 1/2 results in an inequality 

0.01 



4 < 1 (30) 
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that gives a working region in the parameter space of the interferometer. The inequahty 
fl30|) can be expressed in terms of the dimensional experimental parameters as follows: 



where N is the number of atoms in the trap of axial angular frequency 00 and the transverse 
frequency uj±. This inequality gives a fundamental limit on performance of a guided BEC- 
based free oscillation interferometer in Michelson-type geometry. The second fundamental 
limit is due to phase diffusion, which can not be described in the mean field approach. In ad- 
dition, there can be technical limitations like the noise (caused by vibrations), misalignment 
of the splitting laser pulse, etc. 

Figure. O shows the working region of a free oscillation interferometer for a transverse 
trapping frequency u± = 27i x 80 Hz. In the region below the boundary line (which has 
been obtained taking the equality sign in Eq. ( 13T|) ) the interferometric contrast exceeds 50%. 
The maximum number of atoms corresponding to the boundary region for a given trap can 
be read directly from the graph. For example, for u = 27r x 4.1 Hz, N ^ 10^ and the 
interferometric cycle time is 244 ms (the trap period). 

V. DISCUSSION OF SINGLE - AND DOUBLE REFLECTION GEOMETRIES 

In this section, we briefly re-derive results of the analysis of the single - and double 
reflection atom Michelson interferometers [22, [23] and compare their performance with that 
of the free oscillation interferometer. In a single reflection atom Michelson interferometer, a 
BEC sitting at the bottom of a weakly-confining harmonic trap is split into two harmonics, 
which move in opposite directions with the velocities ±Vo. At time r = T/2 where T is 
the interferometric cycle time, a reflection pulse is applied, which adds velocities of =f2Vo to 
the harmonics so that the 1^+ harmonic now moves with a velocity {V+ — 2Vo) and the ip- 
harmonic moves with a velocity (K_ + 2Vo), where V± = ±V are the velocities of the two 
harmonics right before the first reflection pulse is applied. A recombination pulse, identical 
to the splitting pulse is used to recombine the two harmonics at time r = T. 

The total energy of the system is given by Eq. f|T9l) . Unlike in the case of a free oscillation 
interferometer, the energy is not conserved for a complete cycle because of the reflection 
pulses , but remains constant between them. In this case, the velocity of the harmonics at 




(31) 
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the time of recombination can be calculated from Eq. ([T9|) taking into account the momentum 
kick given to the harmonics by the reflection pulse. The difference between the final speed 
of the harmonic and its initial speed AV^ = (Vr — Vq) in the dimensionless variables turns 
out to be: 

AV « (32) 

In deriving Eq. ( |32l) . Eqs. (fT3l) and IHM have been solved to the lowest order in r yielding 
G ~ —t/2 and i? ~ 1 — because the duration of the interferometric cycle in this case is 
much less than the full period of oscillation in the trap. The dimensional version of Eq. fj32|) 
reads: 

Avo ^ jiouTnf, (33) 

where vo = 2hki/M, Lq is the equilibrium size of the condensate and To is the dimensional 
time for the interferometric cycle. Relation (!33|) coincides with the analogous relation (Eq. 
24) in Ref. 27| obtained by a different technique. 

In a double reflection interferometer, two reflection pulses are applied at times r = T/4 
and T = 3T/4 before the recombination pulse is applied at the end of the interferometric 
cycle time t = T. Using the procedure analogous to that described above, we calculate 
the difference between the velocity of the 1^+ harmonic at recombination (Vr) and initial 
velocity (Vq) to be equal to 

or, in dimensional variables. 

This expression matches with the analogous relation (Eq. 43) in Ref. 27|. 

Equations (132 p and have been obtained in the limit when the clouds have zero spatial 
overlap at the time of application of the reflection pulses. 

To compare the performances of the single-, double-reflection, and free oscillation interfer- 
ometers, it is enough to compare the differences between the velocity of the harmonic at 
recombination (Vr) and its initial velocity (Vq) in the three geometries. For a free oscillation 
interferometer, {Vr — Vq) given by the Eq. ( 1271) becomes 
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in dimensional variables. A comparison of the Eqs. (133|) . (135|) and ( 136|) shows that the 
difference in velocities at recombination and initial velocity of a cloud is much larger in a 
single reflection interferometer, less in double reflection interferometer and much smaller in 
a free oscillation interferometer for a given trap frequency. Since this velocity difference is 
the main cause of the loss of contrast (cf. Eqs. (122]) and fl2^ ). an increasingly improved 
contrast can be obtained in a double reflection and free oscillation interferometers compared 
to a single reflection interferometer. 



VI. CONCLUSIONS 



In this paper, we analyzed the operation of a BEC-based free oscillation interferometer 
with optical splitting and recombination of the BEC clouds. Our one-dimensional (ID) 
analytical model is based on the mean field approximation in the Thomas-fermi limit. From 
the ID Gross-Pitaevskii equation, we derive a closed set of ordinary differential equations for 
the parameters describing the shape of the density envelope and the spatially-varying phase 
of the BEC clouds. The derivation is based on the equations of motion for the quantum- 
mechanical expectation values associated with these parameters. The main result of the 
paper is Eq. ( 13T|) . which gives the working region of the interferometer in the parameter 
space and shows how the performance of the interferometer depends on different parameters 
of the experiment such as the number of particles, longitudinal and transverse frequencies of 
the trap, and the velocity imparted by the splitting laser pulses. According to our analysis, 
the reason for the loss of the coherence in a free oscillation interferometer is oscillations of 
the density envelopes of the clouds with a period different from the longitudinal period of 
the trap. 

The analysis of the paper does not include effects beyond the mean field approximation 
such as finite-temperature phase fluctuations along the length of the elongated BEC clouds 



and phase diffusion. Ref. [27| discussed the importance of the phase fluctuations and con- 
cluded that they are negligible for the parameters of the recent experiments [isl, 17, 23 |. 



The phase diffusion, specifically in the context of atom interferometers with the optical 
splitting and recombination of the clouds, has been recently analyzed in j28|. Results of this 
analysis, applied to the case of a free oscillation interferometer, predict that the region of 



15 



good performance is given by the inequality 



(37) 



where uj = {uj^ioY^'^ and a = ^Jh/ (Mu). 



The model of Ref. 



28| goes beyond the mean field approximation by accounting for the 



mode-entangled nature of the two BEC clouds after the splitting, but does not account 
for the development of spatially- varying phases caused by atom-atom interaction during the 
propagation, as opposed to the present paper. Thus, the physics behind Eqs. fl^ and fl57|) is 
complementary, and both these inequalities have to be evaluated and their values compared 
for any particular experiment. 

The relative importance of the effects due to spatially-varying phases caused by atom- 
atom interactions and the phase diffusion is given by the left-hand sides of Eqs. (1311) and 
fl57|) . respectively. The left-hand-side of Eq. ( 13 ip for the parameters of the experiments by 
Burke et al. [15] and Horikoshi et al. [16)] is much less than one, and equals about 0.8 in 
the experiments by Segal et al. 23| . The left-hand-side of Eq. (l37|l for Ref. 



15| is small as 



compared to one and equals to about 0.65 and 1.0 for Refs. [16| and 23|, respectively. This 
shows that the phase diffusion could be partially responsible for the degradation of coherence 
in 16| and that both the effects discussed in this paper and the phase diffusion could be at 



least partial 



Refs. 



15 



16 



y responsible for the loss of contrast in the experiments [23|. The authors of 



23| also list vibrations as a cause for the degradation of the coherence. 



In the experiments discussed in the paper, the frequency of the trap along the guiding 
direction is much less than those along the transverse directions. The BEC clouds are cigar- 
shaped with the largest dimension along the weak guiding direction of the trap and are 
moving along the same direction. This is the reason why a ID theory is a good approximation 
to the experimental situation. A possible slight misahgnment of the optical splitting pulses 
can result in a more complicated 2D or 3D motion of the BEC clouds and their rotations. 
Analysis of such dynamics requires generalization of the ID model of the present paper to 
higher dimensions. 
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FIG. 1: (Color online) A schematic of a free oscillation interferometer. 
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FIG. 2: Time evolution of the dimensionless radius i? of a harmonic and the quadratic phase G 
(rad) for a trap period. The horizontal axis is dimensionless time from to 2it. The inset picture 
shows the effect of interatomic interactions on G when the two harmonics pass through each other. 
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FIG. 3: Time evolution of the cubic phase S'(rad). The horizontal axis is dimensionless time, r 
from to 27r. The cubic phase develops only when the clouds overlap during splitting, when they 
pass through each other and when they recombine. The inset picture shows the evolution of S 
when the clouds pass through each other. 
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FIG. 4: The probability density |'0oP versus the coordinate ^ for A.K = 5 and F = 0. The 
probabihty density oscillates several times under its envelope that reduces the contrast of the 
interference fringes. 
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FIG. 5: Working region in parameter space of a free oscillation interferometer , with the longitudinal 
trap frequency uj (rad/s) and the total number of atoms N in the condensate. The interferometer 
works in the unshaded region and does not work in the shaded region. 
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